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ABSTRACT 


Data  taken  from  the  second  generation  of  the  Naval 
Postgraduate  School's  ocean  floor  aii  land  geomagnetic  data 
collection  system  is  analysed  to  determine  the  wave  nature 
of  Geomagnetic  Hicropulsations  in  tia  .05  *o  13  Hz  range. 
The  measurements  were  naia  in  Honterey,  California  luring 
June  to  September  1982.  Thraa  of  thass  experiments  war? 
twenty- four  hours  in  length.  During  two  such  axoeri^-"*T 
land  and  sea  data  ware  t~ans:nitcei  :o  a  central  point  where 
the  signals  were  racordai  Bimultaneoisly.  3eoaa gnat ic  fluc- 
tuations were  analyzed  to  detar:nin  =  various  factors  charac- 
terizing the  wav=s.  These  factors  include  the  power  spectral 
densities,  coherences,  Stores  para^i: =rs ,  orientation  angle, 
degree  cf  polarization  anl  ellipirit/.  The  software  devei- 
oped  in  this  ra search  oro/ides  an  extensive  capability  for 
analysis  of  geomagnetic  noise.  Long  term  lata  collection  cat 
utilize  it  :o  establish  traits  in  both  tha  spatial  ani 
teaporal  prooerties   of  naturaLly   ooouring  EL? 
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i.   introduction 

This  thesis  is  part  of  the  long  terra  effort  of  the  Naval 
Postgraduate  School  to  callect  and  interpret  ELF 
elctromacnetic  noise  on  tie  sea  floor  and  from  a  nearby  land 
station.  The  project's  objective  is  to  study  md  interpret 
signals  in  the  0.05  to  10  3  z  frequency  ranne.  It  emphasizes 
the  importance  of  obtiiiing  measure  a ents  of  geomagnetic 
noise  en  the  s?a  fleer  anl  on  land  dvbc  a  periDl  of  several 
years. 

The  specific  sbective-  of  this  :hasis  are  twofold;  to 
calibrate,  install,  anl  ccecata  1 1  s  land  data  collection 
system;   and  to   ievelope  a  computer  program  zo    analyze  the 
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A. 


IIHE    VARIATIONS     OF    THE    GEOMAGNETIC    7IELD 


The  earth's  magnetic  field  has  been  used  for  centuries 
for  navigation  at  sea  a.ii  for  sar/ey  lirection  on  lata. 
These  cemmon  ises  are  Impendent  Da  the  static  magnetic 
field.  The  static  field  ran  be  generally  characterized  as  a 
magnetic  dipole  at  a  11.5  iegrss  inclination  id  the  axis  of 
rotation  of  the  earth.  It  shows  ::i/  minor  diurnal  varia- 
tions in  direction.  Its  aagnitude  /iries  greatly  over  the 
earths  surface  but  usually  is  in  the  range  of  43,300  to 
70,000  nano-Teslas  (nT) .  It  Honterey,  California  the  static 
field    has   a    value  of    48,633    nT. 

It  has  cniy  been  siacs  the  discsvsry  of  magnetic  induc- 
tion that  the  underlying  /ariations  La  the  static  field  have 
become      well     known.  rhese   /aria:L:ns      are      flue      to      many 

soirees.      The   sources   iaclide,    but   ar=    not    limited    vo,    solan 


activities,  Lunar  daily  variations,  magnetic  storms,  thun- 
derstorms in  the  earth's  atmosphere,  ear-  hquafces,  and  the 
detonation  of  nuclear  dsvioes.  3ns  o!  the  major  contributors 
to  this  is  the  interactioi  of  the  charged  particles  from  -he 
su.i  with  the  earth's  magna  t  osphere.  rhese  events  vary  widely 
in  duration,  amplitude,  ail  the  freqiency  range  they  effect. 
Within  the  frequency  rang?  of  this  system  the  most  common 
variations   are   characterized   as    uiior d d uisat ions. 

Gecmacnetic  micr  cpuisa  t  ior.s  are  defined  to  be  fluctua- 
tions in  the  earth's  magnetic  field  which  have  periods  f:on 
0.2  seconds  to  1 0  minutes.  They  vary  La  amplitude  from  0.001 
nT  to  as  high  as  a  few  teis  of  nl.  Iney  are  attributed  to 
the  interaction  of  groaps  of  ohargeii  particles  propagating 
in  the  earth's  magnetospher e  as  hydra  magnetic  waves.  These 
waves  ir.turn  create  the  soiree  currents,  (primarily  in  the 
ioasphere),  for  the  mioropulsatio^  field.  Kioropulsatiors 
are  divided  into  two  subcatagories:  irregular  and  conti- 
nuous. The  abr  evi  at  io.i  s  for  these  ar*=  pc  for  pulsation 
continuous  and  pi  for  pdlsation  irregular.  The  primary 
difference  between  the  t*o  is  that  the  pc  mioropulsat ions 
have  a  semi- sinusoidal  nature  while  tie  pi  miooopaisations 
show  irregularities  in  both  amplitude  and  frequency.  The 
period  ranges  and  average  amplitude  for  these  pulsations  ars 
listed  in   Table    I,    [ Ref .    1]. 

The  overwhelming  majority  of  the  pulsations  in  this 
frequency  range  are  due  to  solar  activity.  En  addition, 
"Schumann"  resonances,  will  also  be  apparent  in  the  vicinity 
of  6-8  Hz.  These  resonances  are  due  to  lightening  transients 
in  the  concentric  spherical  cavity  between  the  earth's 
surface   and   the     lower    region   of    the    ionosphere,    *Ref.    2]. 


r ABLE    I 


Geomagnetic:   aicropulsatlon    Clc.sses 


Notation 
PC  2 

P34 

oil 
pi  2 


Period  (Szz ) 

0.  2-5 

5-10 

13-45 

45-150 
150- 600 
1-40 
40-150 


\3tpli4:ude(nr) 

0. 05-0.1 

3.  1-1 

0.1-1 

0.  1-1 

1-10 
0.01-0.  1 
1-5 


SYSTEM    DESCRIPTION 


Figure    1.1        Sensing   Coil    Dimensions 
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Ths  ser.sor  system  c.irren-'riy  employed  consists  of  a  three 
ceil  land  system  deploys!  at  La  Mesa  Village  housing  area, 
approximately  2  km  from  tie  Naval  Postgraduate  School  data 
collection  center,  and  a  t«ro  coil  sea  system  deployed  on  the 
bottom  cf  Monterey  Bay.  Its  distance  to  the  Naval 
Postgraduate  School  (NPS>  varisi  from  2  to  U  k.m.  The 
original  intention  of  the  project  rfas  to  deploy  the  land 
system  approimately  22  nLlss  from  the  Naval  Postgraduate 
Scnool  at  Chew's  Ridge.  Pais  provsl  impossible  due  to  data 
transmission  dif  ficultiss.    Sss  Appendix   3  for   additional 


Each  sensing  coil  is  i  coatia.i3.isly  wound  coil  antenna 
made  up  cf  approximately  54  60  turns  of  18  gaugs  copper  wire. 
Ths  averaqe  sensing  arsa  of  each  coil  is  3.3825  square 
meters  and  its  dimensions  ars  as  ispicted  in  Figure  1.1. 
Ths  sea  system  has  its  coils  mounts!  orthogonally,  usually 
in  the  horizontal  plans.  The  land  system  coils  ars  affixed 
to  a  wooden  mount  that  hoi  is  them  ii  an  orthogonal  configu- 
ration; see  Figure  1.2.  rhe  iirsctional  orientation  of  oh? 
sea  system  cannot  be  determined  when  deployed.  The  land 
system  is  oriented  wita  the  X-coii  in  the  magnetic  north- 
south  direction  ,  y-coil  in  the  magnetic  east-west  dirsc-icn 
and  the  Z-coil  in  the  vertical. 

The  two  systems  were  designs!  to  be  as  similar  as 
possible.  The  coils,  prei mpiif ier,  signal  coiditioner  ani 
pulse  code  modulator  (PCM)  ars  identical  in  design  ?nd  were 
all  calibrated  to  insure  tie  production  of  valid  data.  The 
technical  description  of  aLl  of  these  items  ma/  be  found  in 
[Ref.  3]  .  The  land  system  is  ispLoyed  as  shown  in  Figure 
1.3.  The  pre- amplifiers  are  battsry  powers!  and  located 
within  one  foot  of  the  coils.  Ths  roil  mcuno  and  preampli- 
fiers were  positioned  approximately  50  metsrs  from  tae 
signal  conditioners,   and  connected  by  three  coaxial  cables. 
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Figure    1.2        Sensor    Mounting   Blocie 

Ths  signal  conditioners,  PCM  boari  aad  transmitter  are  all 
housed  in  a  small  wooden.  aiilding  ani  powered  by  D.3.  power 
supplies.  This  building  is  the  closes-  source  of.  alternat .ing 
current. 

After  the  data  has  besn  encode!  ay  the  PCM,  it  is  tran- 
smitted by  7HF  radio  data  link  to  tie  Naval  Postgraduate 
School    for   recording. 

C.       DATA    COLLECTION    AND    PROCESSING 

Land  and  sea  data  are  simultaneously  received  ani 
recorded  on  separate  channels  of  an  analog  tape  recorder. 
Each    tape    holds    60    to    120    ninutes    3f    lata.      These    tapes    must 


12 


, — r 


!        i    V 


•<::  wplifie!  ; 


\TTERY 


PRECCNP1  riON&R 


rRASSMITTER 


'■H 


t 


Figure  1.3    Laai  Data  Collection  Systea 

be  decoded  and  re- recce  del  onto  ligltal  tap?  euitaole  for 
th  =  main- frame  computer  a:  the  Churci  Computer  Center.  Each 
digital  tape  CDn tains  93  nlnutes  of  Land  cr  sea  data.  Zach 
of  these  digital  rapes  is  then  processed  through  th9 
ccaputer  program,  (See  Ciaptec  3f  Section  C),  to  prcviie 
the  plots  of  the  various  interpretive  factors.  See  Reference 


3  for   the  details 


,-p   - 


he      recording,      decoding      and   re-re- 


cocding    en    diaitai    tape. 
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II.   SYSTEM  CALIBRATION 


A.   INTRODUCTION 

The  triaxial  ceil  najaetomstar  system  is  diagramed  in 
Figure  2.1.  The  data  sensed  by  this  aagnetomet  ;r  system  is 
giver,  by  Vx  (t)  ,  Vy(t)  and  7Z  (t|  ;  ths  output  voltages  ger.?r- 
atsd  by  the  input  magnetic  fields.  Tia  system  transfer  func- 
tion must  be  dsterminad  so  that  the  magnetic  field 
conponer.ts  can  be  removed  fron  tha  output  voltage.  In  the 
tine  domain  the  relationship  between  the  output  voltage  and 
tha  magnetic  field  is  a  convolution  Integral, 


v;  (t)=J  a;  (t-f)  a;  (f»  it' 


(D 


whara    i=x,y,z.    It  is    ao:?   sonvanient    to    express    Egn.       (1)     in 
the    frequency    domain. 


Vj    (ui)  =H;  (u^S,-  M 


(2) 


Tha  transfer  function  HI,'  (t*>)  can  easily  be  determined  in 
terms  of  the  ratio  of  output  to  an  Input  test  signal  or  by 
analysis    of    the    measur=neic    system's    properties. 

The  next  section  discisses  tha  ronpenents  of  tie  systam 
and  attempts  to  develop  i  systen  nodal.  The  last  section 
describes  tha  calibration  n easurament s  and  the  experimental 
determination    of   the    transfer   function. 

3.       SYSTEM    DESCRIPTION 

The    individual      transfer    fmctioa      for   the      sensing   coil 


Vc  (cu)  =  jw/K3  («*») 


(31 
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Figure  2.1    Systea  Coaponsnt  Configuration 

where  Vc  {(At)      is  the  spertnm   of  tii  voltage  directly  out  of 
the  coil  and  K  is  i  gain  factor. 

The  r.ex*  component  In  :  he  systsa  is  the  preamplifier  low 
pass  filter.  The  preampLifier  in trD duces  a  gain  of  1000  so 
the  voltage  signal  can  be  sent  to  remote  instrumentation. 
The  low  pass  filter  is  i  five  p:)Lr  Chebyshev  with  a  3dB 
cutoff  at  20  Hz.  The  filter  is  oriiiirily  usei  to  remove  S3 
Hz  signals  frou  the  systen  and  to  pref liter  the  data  before 
it  is  digitized.  An  additional  pole  is  introduced  because  of 
a  oapacitive  coupling  between  the  second  and  third  stages  of 
the  filter.  Figure  2.2  shows  a  pr>l  =  -zero  diagraa  of  the 
pr= amplifier.  The  transfer  function  takes  the  fnum: 

H(j«w)  =K/(jw+a)  (j^  +  b)  (j«u*c-jd)  (jw  +  ovjii  (jo+e-jf>  (ju/+e  +  jf)  ((*} 
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whare  K  is  -he  gain  cf  tis  preamplifier  and  the  coeficisnts 
a  -  f  are  determined  from  the  pols-zsro  diagraa. 

The  prscondi tioner  na  =  additional  gain  of  a pproximately 
1030  and  has  zen  =r  diodss  to  1  i en  1 1.  L.ir  input  voltage  to  the 
PCS  board.  Typically  the  oitoff  is  a::and  7  to  3  volts. 

The  magnitude  and  phise  rasponsa  of  the  oreaaipiifier- 
filter  and  preconditions r  are  shown  In  Figures  2.3  through 
2.3.  The  input  test  signal  generating  these  curves  had  flat 
magnitude  and  phase  character  is  snowr.  in  Figures  2.9.  rhs 
figures  then  represent  tha  actnal  transfer  functions  of  tha 
preamplifier- filter  prsconditioner  system  for  the  x-y-z 
channels  of  land  and  oceai  sensor  systems.  The  figures  show 
that  the  channels  are  well  matched  in   magnitude  and  phase. 


phase  response 


is  iLiear 


an  -^ 


virtually  identical   in 


every  case.  This  is  an  iaportant  point  since  the  phase  paru 
cf  the  transfer  function  is  not  r  shoved  frnn  the  out  out 
voltage  signal.  However,  since  the  lata  analysis  is  in  terns 
of  power  spectral  densities  and  coherences,  this  fact  has  ?.z 
bearing  en  the  results. 

The  total  transfer  function  of  tns  system  is  the  produot 
of  Zgn.  (3)  and  (U)  .  Basel  on  these  formulas  and  tha  figures 
tha  resulting  transfer  fnnction  is  predicted  to  be  fairly 
linear  at  the  low  frequency  end. 


C.   DETERMINATION  0?  TH2  31 STZX    TRANSFER  FUNCTION 

The  magnitude  of  the  system  transfer  function  was  deter- 
mined experimentally  through  tie  use  of  a  HelmhoLtz  coil,  a 
source  of  sinusoidal  cirrant  and  a  spectrum  i.naiyzer.  The 
experimental  set  up  is  illustrated  in  Figure  2.10.  A  one 
turn  Helmholtz  coil  with  a  0.5  1  iirter  radius  and  a  0.51 
meter  seperation  between  ooiis  was  ised  to  provide  a  cali- 
brated magnetic  field  to  tae  sensing  toil.   The  sensing  coil 
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was  placed  concentric  to  the  HslmhoLtz  coil,  and  centered 
vertically,  to  ensure  a  uniform  magnetic  field  about  the 
sensing  coil.  The  magnetic  field  in  this  rsgior.  of  the 
Helmholtz  coil  is  directai  along  its  axis  and  is  given  by 
[Rsf.    41    : 

B    =>^8  N 1/5  T^R  (5) 

to  terms  on  the  crier  of  (L/R),  whara  L  is  the  linear  dimen- 
sion of  the  region,  N  is  the  number  of  turns,  R  is  the 
Helmholtz  ceil  radius  ani  I  is  the  driving  current  of  the 
coil.  The  driving  current  was  produced  by  a  Wavetek  signal 
generator,  Model  Number  142.  The  current  was  measured  by 
reading  the  voltage  across  a  996  ohm  resistor  in.  series  with 
the  Helmholtz  coil.  A  spa  strum  analyzer  (HP-3592A)  was  used 
to  measure  the  system  h:puj-  to  avoid  any  spurious  signals 
generated  by  harmonics.  By  repealing  the  measurements  at  a 
number  cf  different  frequencies  and  drive  currents,  the 
system  transfer    function    oould   tie   measured. 

The  actual  caiibra~iDP.  was  conducted  at  the  La  M.esa 
village  site  using  the  setup  illustrated  in  Figure  2.13. 
The  Helmholtz  ceil,  sensing  coil  and  preamplifier  wars 
located  approximately  40  maters  from  the  shack  which  housed 
the  remainder  of  the  experiment.  Jith  the  Helmholtz  coil 
disconnected  the  background  magnetic  noise  of  this  location 
was  repeatedly  examined  to  ensure  accurate  readings.  Figure 
2.11  shows  a  typical  background  noisa  plot  froa  the  spectrum 
analyzer  (the  transfer  function  has  not  been  rsmovsd)  .  The 
distinct  features  are  Schumann  resonances,  otherwise  the 
background  was  low  enough  and  stabla  enough  to  adequately 
calibra-e  the  system.  Dnly  at  the  Schumann  frequencies  was 
the  background  subtracted  from  tha  nsasured  output  voltage. 
Calibration  runs  were  conducted  from  3.05  Hz  to  20.0  Hz  with 
applied      calibrated    magnetic      field    of      0.02,         0.2    and      1.0 
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Figure    2.10        Systsm   Calibration    Experiment 

Figures  2.12  through  2.16  show  the  resulting 
calibration  carves  of  the  land  and  ocean  magnetometer 
systems.  The  true  transfer  function  is  represented  by  the 
0.32  nT  field.  The  larger  applied  fields  show  saturation  of 
ths  system  at  the  protective  zener  iiodes  in  the  precor.ii- 
tior.er.  At  the  low  freguenoy  end  all  the  systems  are  closely 
matched.  At  1  Hz  an  output  voltage  of  1  v  indicates  a 
measured  magnetic  field  of  0.4  nT.  A  1  nT  field  is  saturated 
above    2    Hz. 

The  computation  of  the  equations,  representing  the 
magnitude  of  the  transfer  functions,  requires  fitting  the 
data  to  the  theoretical  transfer  fmctior.  of  the  previous 
section.         This    would   be    lesirabis    it    the    phase    character    of 
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III.       DATA    EVALUATION 


POWER    SPECTRAL    DENSITY     AND    ALIASING 


The  first  step  in  the  analysis  ot  the  data  collected  was 
ths  determination  of  the  power  speotral  densities  (PSD's) 
for  each  channel  of  the  land  ar.d  sea  lata.  This  was  accom- 
plished through  the  use  of  the  computer  program  developed  by 
J.  ;i.  Schweiger  [Ref.  5].  After  elisor  modi  f  icat  ion  of  this 
program,  approximately  twanty  hours  of  real  tine  data  was 
evaluated.  Samples  of  the  PSD's  are  presented  ia  Figures  3.1 
to    3.5. 

The  overall  value  cf  the  PSD '3  is  in  agreement  with 
previous  work  conducted  in  this  area.  The  PSD's  have  a 
decline  of  appr cximately  -6  d3  per  octave  from  0  to  5  Hz 
where  they  tend  to  decrease  their  rate  of  decline.  Night 
PSD's  tended  to  be  3  to  7  is  below  tie  daytime  level.  This 
coincides  with  the  findings  of  McDevifct  and  He nan  (1930),  M. 
fif.  Beard  (1981)  and  A.  C.  Fraser-3mith  and  J.  L.  Buxton 
(1975),    [Ref.    6]. 

One  feature  of  the  PSD's  that  liffered  from  previous 
PSD's  is  the  spike  at  4  HZ  in  the  Laia*  data.  This  was  found 
to  be  "aliasing"  of  the  53  Hz  power  grid  component  by  the 
PCS  sampler.  This  experiment  was  oonducted  with  a  sample 
rate  of  32  Hz;  which,  by  the  sampling  theorem,  should  give  a 
clear  view  of  the  frequency  range  of  0  to  10  Hz,  when 
coupled  with  the  20  Hz  cut-off  filter  in  the  amplifier ,[ Ref . 
7  ].  The    60      Hz      component   of      the      geomagnetic    field      was 

measured  by  ?.  W.  Clayton  (1979)  and  found  to  be  U  (r.T)  /Hz, 
which  is  10  higher  than  the  surrouniir.g  data  "Ref.  8].  The 
20  Hz  filter,  built  into  the  amplification  system,  would 
normally    eliminate    freguenoies   in    this      range    bat    because    of 
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Figure    3.2        PSD,    Y-Coil,    Land   Data 
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the  great  differences  in  magnitude  this  does  nor.  happen. 
When  the  32  to  64  Hz  range  is  convolved  onto  the  0  to  32  Hz 
range  the  60  Hz  component  will  appear  at  4  Hz.  This  destroys 
the  usefulness  of  the  data  in  a  small  region  about  this 
poi  nt . 

Another  point  of  interest  is  illustrated  by  the  vertical 
component  of  the  sea  data  recorded  on  September  16,  1982, 
Figure  3.5.  During  this  experiment  the  Y-coii  was  oriented 
in  a  vertical  rather  than  the  more  normal  horizontal  direc- 
tior..  Examination  of  the  data  shows  13  to  15  d3  decrease  in 
the  PSD.  This  is  in  agreement  with  the  theory  that,  because 
of  the  ocean-air  interface,  geomagnetic  waves  nropagate 
vertically    in    the  sea,    this    eliminating    this    component. 

3.       STOKES     PARAMETERS    AND    THEIR    INTERPRETATION 

After  the  collection  of  the  data  during  Jane  to 
Seotember  of  1982,  the  problem  of  ho*  to  interpret  if  became 
a  prime  concern.  During  the  latter  part  of  this  period  a 
computer  program  was  developed  to  determine  the  Power 
Spectral  Densities  (PSD's)  by  J.;-!.  Schweiger,  Dctober  1982, 
[Ref.    5].      ' 

This  program  provided  the  basic  structure  to  read  the 
data  points  off  of  magnetio  data  tape  and,  through  the  use 
of  a  Fast  Fourier  Transforn  (FFT)  roitine  and  transfer  func- 
tions, to  compute  the  geomagnetio  fluctuations  in  the 
frequency    domain.  This    program    provided   the      basic    format 

which  was  expanded  and  modified  to  produce  other  parameters. 
The  derivations  of  the  equations  for  these  parameters  and 
their    interpretation    comprise   the    renainder   of   this    section. 

The  first  parameter  to  be  discussed  is  the  coherence 
between  the  different  ohanneis  on  land  and  sea.  Coherence 
gives  a  measure  of  the  "singleness"  of  the  source  at  the 
various    frequencies.      "Singleness"    in    this    context    refers    to 
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what  percentage  cf  the  intensity  at  a  frequency  can  be  said 
to  come  from  a  single  source.  To  coapute  the  ooharer.ce  one 
first  calculates  the  cross  power  spectral  densities  between 
the  channels.  For  a  given  pair  of  channels  X  and  Y,  the 
cross  power  spectral  density (S  )  is  iafined: 

N 

SJky(fj)=1/N^X(fj)Mfj     ) 

n=1 


In    -his       aquation  N       is    ths       r.umbsr    Df       blocks    d: 

through   -he   FFT    snd    averagsd;      X(rj)     and    l(fj)    are   the   FFT's 

along    -he    appropriate   axis.         The      zz  hsrence  (£?>       is    defined 

as: 


z 


SyrC^j  ) 


(f;)=- 


*r  v 


iSx  {lj  )  5rifj  )  ) 


For  the  land  data  coherences  wars  also  determined  for  the 
Y-Z  plane  and  the  Z-X  plans.  Figures  3.7  and  3.3  show 
typical  land  and  sea  coherences,  with  all  coils  oriented  in 
tha  horizontal  Diane.  NDte  that  on  the  land  the  coherence 
stays  very  close  tc  its  maximum  vaiua  of  one  until  approxi- 
mately 7-8  Hz  and  then  starts  to  decrease  while  in  the  casa 
of  the  sea  data  there  is  ncra  structure  at  the  lower 
frequencies.  This  type  of  behavior  is  typical  of  most  of  tha 
data,    although   the    exact    frequancy    ol    the    decliia    varies. 

In  addition  to  coherence,  Stokes  parameters  may  be 
computed  from  the  au tc-spactrai  dansitias  and  cross-spectral 
densities.  They  can  be  used  to  interpret  tha  geomagnetic 
waves  on  a  bi-Linear  basis,  a  bi-CircuIar  basis  or  an  ellip- 
tical   basis. 

Although  other  notations  are  utilized  for  the  Stokas 
parameters,  the  most  common  notation;  sO ,  s1,  s2,  s3;  will 
be  employed  here  [Ref.  9],  Tha  general  form  of  the  Stokes 
parameters    are: 
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sO=<a12>+<a2Z> 
s1  =  <d12>-<i22> 
s2=2<a1  a2   :os^  > 
s3=2<a1 a2   sin£  > 


whsre  a1  and  a2  are  tha  instantaneous  amplitude  of  tha 
orthogonal  components  of  the  magnatic  vector  an d  6  is  their 
phase  difference.  In  tsris  of  the  auto-  and  tha  cross-spsc- 
tral    irr.siti?s   the    Stckas   parameters    can    be    written: 

s0=2  [Sx  +Sy) 
s1  =  2(3Jr-S/» 
s2=4  3s(S,y) 

It  must  be  noted  that  tha  Stokes  parameters  ware  developed 
originally  to  describe  monochromatic  light  and  can  only  be 
applied  to  cur  problem  after  the  lata  has  beer,  transformed 
to    the    frequency    domain. 

The  bi-Linear  basis  ir.t  erpretation  of  the  Szokas  parame- 
ters can  be  seen  by  examining  a  superposition  of  two  orthog- 
onal linearly  polarized  waves  [Raf.  10]-  Computing  *:he 
Stokes    parameters  in    this    basis    yields: 

s0=1/2  (B**  Ey  ) 
s 1=1/2  (B^-bJ) 

s2=    3X  By  CDS  {QK  -9y  ) 
s3=    Bx  By  sin  (Qj,  -9X  ) 

Another  basis  usafal  for  interpretation  is  the 
bi-Circular  basis.  This  basis  is  represented  by  the  superpo- 
sition of  two  counter-rotating  circularly  pclarizsd  waves. 
In   this    basis   the  Stokes    para  maters    ara: 


s0  =  B+<-3_ 
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Figure  3.9       Basis   for   Interpretation 

s1=23^.  B_cos  (^  -j>_  ) 
s2=23^B_sin  (£+-£.  ) 
s3  =  3*-B^ 

A  final  basis  Df  interpretation  is  the  elliptical  basis. 
In  this  case  the  basis  is  modeiei  by  the  major  and  minor 
axis  vectors  and  the  angla  formed  between  the  X-axis  and  the 
major  axis  of  the  ellipse.  From  this  basis  the  Stokes  param- 
eters   are: 

s0=1/2  (B^B-J 

s1=1/2  (Bm-EU)  00=    2  <*. 
s2=1/2  (3m-b„)  sin    2  *c 

s3  =  3^  3^ 
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The  important  point  to  note  fron  these  different  basis 
sets  is  that  for  a  given  wave  form  (linear,  elliptical,  or 
cicular)  the  Stckes  parameters,  as  computed  from  the  spec- 
tral densities,  are  the  same  in  any  basis  in  which  they  are 
computed.  This  can  be  easily  visualized  by  taking  a  circu- 
larly cr  linearly  polarized  wave  and  determining  the  value 
of  the  Stokes  parameters  i.i  each  basis.  It  seems  that  these 
interpretations  are  most  easily  visualized  in  the  elliptical 
basis.  Consequently,  in  addition,  to  the  S-.okes  parameters 
-he  ancle  <x.  was  compute!  to  make  the  interpretation  easier. 
This  angle  was  calculated  from  the  Stokes  parameters  s1  and 
s2    by    the   equation: 

©<=1/2    arctan  (s2/s1) 

Examples  of  the  results  of  the  calculations  of  the  stokes 
parameters  and  the  angle  &C  for  land  and  sea  data  are 
displayed  in  Figures  3.10  to  3.19.  Sfote  that  s1,  s2,  and  s3 
are  normalized  by  dividing  by  sO .  In  addition,  for  the  land 
data,  the  same  graphs  were  produce}  for  the  ?-Z  plane  and 
the    Z-X    plane. 

The  use.  of  Stokes  parameters  is  based  en  monochromatic 
and  polarized  waves.  To  ietermine  the  applicaoility  -o  our 
data  of  the  Stokes  paramsters  the  degree  of  polarization 
must    be    found.    For    a   perfectly    polarized    monochromatic    wave 

.2         .1         J.         .,2. 

sO    =s1     +s2    +s3    . 

If  the  wave  is  partially  polarized  then  sO  becomes  greater 
than  the  sum  of  the  squares  of  the  other  three  [RSF.9].  To 
determine  how  close  to  the  polarized  case  the  lata  lies,  the 
degree  cf   polarization    (P)     can    be    calculated.      It    is    defined 

<3l**s2**s3a)* 

P=    . 

s3 


45 


r 


r-i 

S       *occ 

*sooJ                 "Woo-j           -ae-^oo- 

^OG- 

■st 

■        :        i 

, 

:                      -»           : 

wv        . 

cT 
0 

_^g£~                                  i 

g 

| 

^ 

Q. 

zr 

0 

> 

1       i 

i 

/ 

Pi 

i 

ifi_ ; 

\ 

O 

1 
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Figure    3.11        Land    S1/S0,    X-Y    Plans 
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Figura  3.14        Land  &ngl=©<,    X-Y   Plana 
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Figure    3.15        Sea    SO,    X-Y    Plane 
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Figure    3.17         Sea    S2/S0,    X-Y    Plana 
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Figure    3.18         S~a    S3/S0r    X-Y    Plana 
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If  P  is  cne  then  the  wave  is  perfectly  polarized  and  any 
analysis  based  on  the  Stokes  parameters  will  be  completely 
accurate. 

For  values  of  P  less  than  one  the  analysis  is  corres- 
pondingly less  valid.  However,  factors  that  depend  only  on 
s1,  s2,  and  s3,  such  as  the  angle  *<  and  the  ellipticity,  are 
net  effected  by  low  values  of  ?.  la  such  regions  the  noise 
level  is  much  higher  and  the  data  is  not  easily  interpreted. 

The  land  and  s~a  lata  vary  In  their  values  cf  ? 
depending  en  the  time  of  lay  and  otaer  unknown  factors.  In 
general  they  both  shew  a  high  degree  of  polarization  through 
3  Hz  while  above  this  frequency  a  narked  drop  off  occurs. 
If  was  also  noted  that  in  the  land  data  the  degree  cf  polar- 
ization is  virtually  identical  in  all  three  planes.  This  is 
apparently  due  to  the  fact  that  the  polarization  was  not 
oriented  along  any  cf  the  three  sensing  axis  and  therefore 
the  proportion  cf  the  total  intensity  that  resulted  from 
polarized  waves  remained  relatively  constant.  (See  Figures 
3.2  0-3.23) 

Cne  final  tool  cf  interpretation  of  the  lata  is  the 
eliipticify  of  the  wave  (E)  .  The  eliipticity  is  defined  by 
the  majcr  and  minor  axis  cf  the  ellipse ,[ Ref.  3]. 
Specifically  E  is  defined  by: 

tap.  E=  EW/^m 

This  can  also  be  written  in  terms  oe  the  Stokes  parameters 
as: 

i      2     2     K 

E=V2    arcsin  (s3/(s1    +s2   +s3    )     ) 

where  E  varies  from7T/U  to  -  ^"4 .  Th=  interpretation  of  E  is 
illuminated  by  considering  its  dependence  on  s3  and  the 
elliptical  basis.  If  S  is  zero  thai  the  minor  axis  is  zero 
which    means   the    polarization    of   the    rfave    is    linear,    at    E="tf/4 


57 


Figure  3.20    Degree  of    Polarizatl dq ,  Land,  X-Y  Plane 
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Figure  3.21    Degree  of  Polarization ,    Land,  T-Z  Plane 
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Figure  3.22    Degree  of  Pclarizat i on ,  Land,  Z-X  Plans 
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Figure  3.23    Degre-  of  Polarization ,  Ssa 
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the  polarization  is  circular  and  is  rotating  in  a  right- 
handed  sense,  and  at  E=-"/'4  tha  wave  is  again  circular,  but 
with  a  left-handed  rotation.  In  between  these  values  of  5, 
the)  polarization  of  the  wave  is  an  ellipse  with  the  appro- 
priate handedness  of  the  po larization  depending  on  whether  E 
is  positive  or  negative.  The  smaller  the  absolute  value  of 
the  eliipticity  the  narrower  or  more  linear  the  ellipse 
becomes. 

From  F 1  n  u  r  e  s  3 .  2  ^  and  3.25,  the  form  of  these  waves  for 
both  lard  and  sea  data  is  generally  a  narrow  ellipse  with 
varying  left  and  right  handedness  in  the  areas  of  high 
polarization  while  appearing  to  be  nore  elliptical  outside 
these  areas  [Ref-  9]. 

The  final  result  of  this  investigation  is  that  through 
the  computation  of  the  Stakes  parameters  the  nature  of  the 
geomagnetic  wave  can  be  letermined.  Specifically,  to  find 
the  type  of  polarization  or  the  wave,  one  computes  the 
eliipticity  ana  the  orientation  angle  c<  .  rhe  degree  of 
polarization  is  also  needed  to  establish  the  degree  of 
validity  of  the  interpretation.  Naval  Postgraduate  School 
Technical  Reoort  *NPS  61-33-005,  contains  additional  plots 
of  all  the  factors  described  in  this  section. 

C.   COMPUTER  PROGRAM  DE7 ELD PEMENT  AND  DATA  PROBLEMS 

As  noted  earlier  the  ability  to  process  large  amounts  of 
data  in  a  timely  manner  was  first  put  to  use  in  this  project 
in  September  1982.  The  original  program  has  been  expanded 
to  compute  all  of  the  factors  listed  in  the  previous 
section.  It  consists  of  a  main  program,  a  FFT  subroutine, 
and  a  read  subroutine.  &s  noted  before  the  FFI  subroutine 
transfers  the  data  into  the  freguency  domain  and  the  read 
subroutine  reads  the  PCM  lata  from  the  digital  computer  tape 
and  places  it  in  the  input  array.   All  changes  and  additions 
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made  were  in  the  main  program,  therefore  only  the  main 
program  is  included  in  Appendix  A.  This  program,  Landxyz,  is 
usad  to  evaluate  only  t'aa  land  data.  The  version  of  the 
program  used  to  evaluate  the  sea  data  is  identical  to 
Landxyz  with  the  following  exceptions:  sea  system  transfer 
functions  are  used;  bsciuse  tha  s  =  a  system  has  only  two 
coils,    only   one    plane  of    interpretive    factors   is    computed. 

Although   the   basic      program    provad    adequate    to      find  the 
power      spectral    densities,         the   growth      o 


-  v 


the  program 


ue  other  calculations  created  several  problems.  The 


ir.c_u  ee 

Job  Control  Language  (JCL)  orovid=s  the  guidance  to  tha 
operating  system  of  the  computer  to  tell  it  hew  to  run  batch 
jobs.  Specifically  what  input  and  output  devioas  to  use  and 
hew  much  memory  to  allocate  to  tha  program  during  each  step 
of  its  execution.  (In  the  program  tha  JCL  can  be  identified 
by  the  //*s  at  the  start  of  each  lina  cf  this  code.)  During 
the  expansion  of  the  program  the  amounts  of  nemory  in  the 
LIN'K,  GO  and  PLOT  steps  in  the  execution  cf  the  program  have 
all  expanded.  The  currant  JCL  is  caoable  of  handling  up  to 
two    hours    of   data  at    32    Hz    or   one    hour    of    data    at    64    Hz. 

Sea  and  land  data  are  initially  recorded  on  analog  tape 
and  then  decoded  and  re-ra corded  onto  digital  tape  in  the 
Near  Shore  Laboratory  of  tie  Naval  Postgraduate  School.  This 
procedure  is  discussed  in  [Eef,  3].  The  output  of  this 
re-recording  is  approximately  ninty  minutes  of  real  time 
data,  from  either  the  sea  or  land  site,  on  a  digital  tape 
with    a      density    of      800    bits    v= r      inoh     (BPI) .  The    storage 

density  proved  a  great  handicap  in  the  rapid  expansion  of 
the  program.  The  Church  Zomputsr  Canter  has  fifteen  tape 
drives,  but  only  one  of  these  is  capable  of  handling  a  800 
BPI  tape.  The  competition  for  this  one  tape  drive  led  to 
processing  times  of  up  to  twenty-seven  hours.  This  delay 
made    any    programming    error    extremely    costly    in   terms    of    time 
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lost.  An  attempt  was  made  to  transfer  selected  rapes  into 
the  memory  of  the  computer,  but  it  was  unsuccessful  due  to 
minor  recording  errors  an  the  digital  tapes.  These 
recording  errors  were  overcome  by  the  read  subroutine  in 
order  to  be  able  to  work  with  the  data.  Time  did  not  permit 
the  presuit  of  this  avenue  of  Lmproveient  but  an  best  only  a 
few  tapes  could  be  storsd  this  way  because  of  the  vast 
amount  of  raw  data. 

The  computer  proaran  rfas  validated  through  the  use  of  a 
test  tape.  The  test  tape  was  recorded  with  a  sinusoidal 
input  at  known  frequency  and  voltage  level  as  measured  by  a 
soectrum  analyser.  when  this  tape  was  processed  by  the 
Landxyz  program,  the  correct  frequency  and  amplitude  were 
produced,  [P.ef.  11]. 

A  noise  floor  tape  was  produced  by  replacing  the  sensing 
coils  with  120  ohm  resistors.  This  tape  was  then  evaluated 
by  the  Landxyz  program  to  ore  luce  ths  PSD  of  the  noise  fiscr 
of  the  system.  Figure  3.25  shows  this  noise  floor  superim- 
posed en  a  land  PSD  plot. 
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IV.        CONCLUSIONS    AND    RE CD  EMENDATIONS 

A.       CONCLUSIONS 

As  stated  earlier  part  of  the  objective  of  this  thesis 
was  to  calibrate,  install  and  operate  the  land  portion  of 
the  Naval  Postgraduate  School's  gesiagnet ic  data  collection 
sys-em.  This  was  generally  successf  ul  ly  accomplished  with  the 
exception  of  degradation  of  data  produced  by  using  the  more 
noisy  La  tfesa  site  rather  than  the  aore  remote  Chew's  Ridge 
site.  Exactly  how  much  of  an  improvement  in  the  quality  of 
the  data  will  result  from  the  move  to  Chew's  Ridge  is  diffi- 
cult to  predict.  It  could  possibly  diminish  the  50  Hz  compo- 
nent so  that  data  could  be  obtains!  in  the  region  of  u  Hz. 
Overall  some  70  hours  of  land  lata  wis  collects!  during  this 
research,  of  which  appro xi lately  twenty  hours  has  been  eval- 
uated to    date. 


For    the   thesis's    second    objective. 


r  Vie, 


production  of  a 


computer  program  and  procedures  to  interpret  the  data,  the 
expansion  of  the  basic  program  to  oompute  the  coherences, 
Stokes  parameters,  degree  of  polarization,  elliptical  orien- 
taion  angle,  and  the  ellipticity  would  appear  to  have  satis- 
fied the  reguire  ment .  The  major  drawback  to  this  system  is 
the  limited  amount  of  lata  that  can  be  processed  in  a 
reasonable  time.  At  most  the  current  programs  ars  capable  of 
evaluating  only  5-6  hours  of  land  and  sea  data  in  a  twenty- 
four  period.  This  is  primarily  due  to  +ape  drive  avail- 
ability, but  the  requirement  of  the  man-hours  for  decoding, 
re-recording  on  digital  tape  and  filing  these  tapes  at  the 
computer  center  cannot  be  discounted.  The  data  processing 
system  has  reached  the  stage  where  large  amounts  of  data  can 
be  evaluated  but  it  is  not  yet  ready  to  operate  on  a  conti- 
nuous  basis. 
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3.       RECOMMENDATIONS 

The  first  recommendation  is  that  the  data  transmission 
problem  be  solved  so  that  the  Shew's  ridge  sita  can  be  oper- 
ated on  a  continuous  basis.  This  has  two  advantages;  in  the 
improvement  in  the  guality  of  the  sampling  environment  and 
in  giving  a  contiuous  data  stream  that  can  be  recorded  at 
appropriate  intervals  to  provide  a  Long  term  data  base  from 
which    lew    freguency    trends    could    be    observed. 

Improvement  could  also  be  mads  in  the  ansa  of  data 
processing.  Specifically,  a  new  method  of  data  handling  is 
needed.  Several  alternatives  exist,  the  first  being  the 
acquisition  of  the  projects  own  decoder  and  digital  tape 
drive  to  provide  a  digital  tape  of  a  higher  density  than  can 
now  be  made.  Phis  would  T.ake  it  possible  for  a  two  to  eight 
fold  reduction  in  the  number  of  tap^s  necessary  to  hold  the 
data  and  weald  also  giv?  faster  turn — around  time  during 
processing  due  tc  the  increased  number  of  tape  drives  avai- 
lable at  the  computer  center.  Other  data  handling  schemes 
are    also      feasible    if   the    problem      of    how    to    rsad      data    from 
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digital   tapes  into   the   mass   storage  system 


the 
computer  can  be  overcome.  This  would  make  it  possible  tc  run 


programs  in  other  than  the  "batch"  mode  and  avoid  the  use  of 
tape  drives  to  gain  access  tc  the  data. 
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APPENDIX    A 
COMPUTER    PROGRAM 

//LANDXYZ    JOB     ( 1  573 , 0 129)  ,  •     POGOE  SMC    1  573' , SLA SS  =  F 

//♦MAIN    0RG=NPGVM1.  1  573P,LI  NES=(99| 

//♦FORMAT    PR, DDN  ftME^PLOT . 3 Y  SVECTR, DE3  T  =  LOC AL 

//    EXEC    FRTXCLGP,REGION.  GD=  3000 K 

//FORT.  SYSIN    DD     ♦ 

INTEGER*2     IN  (16) 

C  ARRAY    'IN'    IS    USED    IN    READING    DATA    FROM    TAPS 

COH?LEX*8     XX  (3  192)  ,YY  (3192)  ,  Z  (3  19  2) 
COMPLEXES     SXY  (  3  192)  f  S  YZ  (3  1  9  2)  ,  SZX  (8192) 

C  THE    COMPLEXES    ARRAYS     ARE    USED    13    ORDER    INPUT 

C  DATA    AND     INITIALLY    REPRESENTS    VOLTAGE-TIME 

C  SERIES    INFORMATION 

DIMENSION    TIME  (3192)  ,  FREQ  (8192k  ,SORK  (16384)  rFRQ2  (3192) 

DIMENSION     IX  (3  1  92)  rZY  (8  19  2)  ,ZZ(3192) 

DIMENSION     SOXY  (31  9  2)  ,  SO  YZ  (3  1  9  2)  ,  SOZX  (8192) 

DIMENSION     S1XY  (819  2)  ,  S1YZ  (3192|  ,  S1ZX  (8192) 

DIMENSION     S2XY  (3192)  ,  S2  YZ  (8192»  ,S2ZX  (8192) 

DIMENSION     S3XY  (8192)  ,  S3Y2  (8192)  ,  S3ZX  (8192) 

DIMENSION    COXY  (8192)  ,  COYZ  (8192)  ,  COZX  (3192) 

DIMENSION     AXY  (8192)  ,  A  YZ  (3  19  2)  ,  AZX  (819  2) 

DIMENSION     PXY  (8192)  ,  ?  YZ  (8  19  2)  ,  ?ZX  (319  2) 

DIMENSION     DEXY  (8192)  ,  DEYZ  (8192i   ,DSZX  (3192) 

C  THE    'Z'     ARRAYS    REPRESENT     FREQUENCY    ECMAIN     (EFT) 

C  MAGNITUDE    DATA    AND    ARE    EVENTUALLY    CONVERTED    TO 

C  POWER    SPECTRAL    DENSITY    INFORMATION.     COXY    IS    COHERENCE 

C  MEASURED     X   TO    Y.     AXY     IS    THE    ANGLE    BETWEEN    THE    MAJOR 

C  AXIS    OF    THE    ELLIPSE     AND    THE    X    AXIS.     PXY    IS    THE    DEGREE 

C  OF    POLARIZATION    IN    THE    X-Y    PLANE. 

C  DEXY    IS    THE    ELLIPTIIITY    IN    THE    X-Y    PLANE. 
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INTEGER*!     IT 3  (  12)  /12*  0/ 

REAL*U  RT3  (28) /2S*0.0/ 

REAL    ALAB  (  11)  /'X'  ,  '  Y«  ,'Z' ,,C',,30,,,S1,,,32,,'S3',,A' 

r«P<  ,  'EV 

REAL*8   TITLE  (12) 

EQUIVALENCE(TITLE(1>  ,  RT3(5)  ) 
C  ARRAYS     •  ITB»f  »RTB*r  •  ALAB1,  AWD    TITLE1     ARE    USED    IN 

C  THE    VERTS  ATEC    PLOTTER    OUTPUT. 

DATA    XX,YY/16384*(0.  ,  0.  )  / 

DATA    Z/81  9  2*  (0.  ,0.  )  / 

DATA    ZX,ZY/16334*0./ 

DATA    COXY,  SOXY/16334*  0./ 

DATA    52XY, S3XY/163 34  *  0.  / 

DATA    S1XY,  S1YZ/16384*0./ 

DATA    COYZ,  S0YZ/16334*0./ 

DATA    S2YZ,  S3Y Z/1638^* 0. / 

DATA    COZX,  S0ZX/1638U*0./ 

DATA    S2ZXrS3ZX/16384*0./ 

DATA    ZZ,?REQ/16384*a.  / 

DATA    AXYr  AYZ/16384*D.  / 

DATA    AZX, S  1ZX/ 16384*0./ 

DATA    PXY,  PYZ/16384*D.  / 

DATA    PZX,DEXT/16384*0  ./ 

DATA    DEZX, DEYZ/16384*  0.  / 

?I2=6.2831 853 
C  THE    FOLLOWING    SECTION    READS    THE    DECODED    PCM    DATA    OFF 

C  THE    DIGITAL    TAPE    AMD    FOR    EACH    CHANNEL    SKIPS    EVERY 

C  OTHER    SAMPLE    POINT. 

C  I5EC=3 

C  ITL=TSSC*3  2 

ITL=300 

DO    55    JJ=1 ,ITL 

CALL    RD(20,IN,200,IREC,ISR) 
55    CONTINUE 
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IFRAME=8192 

NR=16 
FNR=FLOAT  ( NR) 

DO    70    L1=1  fNR 
C  THE    DO    LOOP    ENDING    9 ITH    STATEMENT    70    ENABLES    THE    PEO- 

C  GRAM    TO    PROCESS    A    LA RGE    AMOUNT    OF    DATA    BY    REPEATING 

C  THE    PROCESS    IN    BLOCKS.    THE    DATA    POINTS    FROM    EACH    RrJN 

C  THROuGH    THE    DO    LOOP     ARE    ADDED    TOGETHER    AND    EVENTUALLY 

C  AVERAGED     EY    THE    THE     NUMBER    3?    TIME    THROUGH    THE    LOOP. 

C  'NR'     REPRESENTS    THE     NUMBER    OF    DATA    SEQUENCES    TO    BE 

C  AVERAGED 

C  1    SEQUENCE   CURRENTLY     EQUALS    3192     DATA    POINTS    FOR    EACH 

C  CHANNEL    OR   256    SECONDS    AT    A    SAMPLE    RATS    OF    32    CPS. 

DC    6  0    JJ=1  fIFRAME 

CALL    RD  (20  ,IN,  1000  rIR  ZCf  IRR) 

XX  (J J)  =117  (2) 

YY  (JJ)=IN  (3) 

Z(JJ)  =IN(4) 
6  0    CONTINUE 

WRITE  (6,200)  IRRrIREC 
2  00    FORMAT ( 10 X,' IEP=' , 15 , 5X, ■ IREC=' , 16 , /) 
C  THIS    SECTION    GENERATES    THE    TIME    AND    FREQUENCY 

C  ARRAYS    AND    NORMALIZES    THE    INPUT    PCM    DATA    TO    VOLTAGE 

C  IN    PREPARATION    FOR     FAST    FOURIER    TRANSFORM    TO    THE 

C  FREQUENCY   DOMAIN. 

N=819  2 

FN=FLOAT(N) 

DELTAT=1./32. 

T=FN*DELTAT 

DELTAF=1.  0/T 

DO    20    J=1 , N 

TIME  (J)  =  D  E  LT  AT  *FL0  AT  (  J) 

FREQ  (J)  =DELTAF*FLOAT(  J) 

XX  (J)  =  (XX  (J)  -2  045.5)  *  10. /2  0U  5.  5 
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XX  (J)  =REAL  (XX  (J)) 
YY  (J)  =  (YY  (J)  -2045.  5)  *  10./2  045.5 
YY  (J)  =REAL  (YY  (J)) 
Z  (J)  =  (Z  (J)  -2  04  5.5)  *10  ./20  4  5.5 
Z  (J)  =  REAL  (Z(J)  ) 
C  'XX'     IS     THE    X-COIL    DATA,     •YY'     15    THE    Y-COIL    DATAr 

20  CONTINUE 

WRITS  (6,9  0  0)  X,  Y,Z? 
201  CONTINUE 

DO  21  J=1 , N 

PBQ2  (J)  =ALOG10  (FSSQ(J)) 

21  CONTINUE 

C        THE  NEXT  3  STATEMENTS  PERFORM  AN  ??T  ON  THE  INPUT 
C        TIMS  SERIES  DATA. 

CALL  FOOST  (XX,  N,  1,-1 ,  0,  WORK) 

CALL  FOUHT  (YY,N,1,-1,  0,HORK) 

CALL    FOORT  (Z,  H,  1,-1,3  ,WORK) 
C  THIS    BLOCK   OF    STATEMENTS     APPLZ     IKS    SYSTEM     (VOLTAGE    13 

C  B- FIELD)      TRANSFER    FJNCTION    ID    THE    TRANSFORMED 

C  FREQUENCY    DOMAIN    DATA.    THIS    BLOCK    ENDS    AT    9. 

C  THE    TRANSFER    FUNCTION    CONVERTS    VOLTS    TO    NANOTESLAS. 

C  *****WARNING***** 

C  THIS    TRANSFER    FUNCTION    IS    GOOD    FOR    MAGNITUDE    ONLY. 

DO    9    L=1fN 

FRQ=FREQ(L) 

IF  (FRQ.LE.  25.)  GO    TO    1 

XX  (L) =XX(L)/29. 

GO    TO    3 

1  IF  (FRQ.  LE.  15.  )  GO    TO    2 

XX(L)  =XX(L)/(1  05.5-3.  14*FRQ) 
YY(L)  =YY  (L)/(181.3  2-7.588*FR2) 
Z  (L)=Z  (L)  /  (17  7.26-7.  '4  84*FRQ) 
GO    TO    8 

2  IF  (FRQ.LE.  10.)  GO    TO    3 
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XX  (L)=XX(L)/(5.958*F3Q-30.97) 

YY  (L)  =YY(L)/(7.  166  *FR  Q-39  .  99} 

Z  (L)  =Z(L)  /  (6.4  9*FRQ-3  2.35) 

GO  TO  8 
3    IF(FRQ.LE.  7.5)  GO    TO    '4 

XX  (L)=XX  (L)/(3.492*FRQ-6.  31) 

YY  (L)  =YY(L)/(4.252*FRQ-10.35| 

Z  (L)  =Z  (L)  /  (4.  044*FRQ-7.89) 

GO  TC  3 
a    IF  (FP.Q.LE.  5.)  GO    TO    5 

XX  (L)=XX  (L)/(2.631  1*PRQ+0- 1466  7) 

YY  (L)  =YY(L)/(3.012*PRQ-1.  5  5) 

Z  (L)  =Z  (L)  /  (3.  184*FRQ-  1.44) 

GO    TO    8 

5  IF  (FP.Q.LE.  3.)  GO    TO    5 

XX (L)  =XX(L)/(2.631 1 *F RQ  +  0 .  145 6 7 ) 

7  YY  (L)=YY  (L)/(2.702*FRQ) 
Z  (L)  =Z  (L)  /  (2.  92*FRQ) 

GO    TO    3 

6  XX  (L) =XX(L)/(2.72*FR2) 
GO    TO    7 

8  CONTINUE 

9  CONTINUE 

C  THIS    BLOCK,     ENDING    WITH    30,    T3E     MAGNITUDE    OF    THE 

C  FREQUENCY    DOMAIN    DAI  A    13     TAKES     AND    PLACED    IN    A    REAL 

C  ARRAY, 'ZX'    FOR    THE    X-COIL,'ZY'     FOR    THE    Y-COIL    AND    ZZ 

C  FOR    THE    Z-COIL.    THIS     DATA    IS    DIVIDED    BY    THE    NUMBER    OF 

C  SAMPLE    POINTS    AND    SQUARED    TO    DETERMINE    POWER. 

DO    30    11=1  ,N 

ZX  (I1)=ZX  (11)  +  (CABS  (XX  (11)  )  /FN)  **2 
ZY  (11)  =ZY  (11)  +  (CABS(YY(I1)  )  /FN)  **2 
ZZ(I1)  =ZZ  (11)  *  (CABS(Z  (11)  )/FM)  *  *  2 
SXY  (I1)=SXY(I1)  +  ((XX  (11)  /FN)  *C0NJG  (YY  (11)  /FN)  ) 
COXY(I1)=COXY  (11)  ♦CABS  (  (XX  (11)  /FN)  *CONJG(YY  (11)  /FN)  ) 
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SYZ  (11)  =SIZ(I1)  ♦  ((YY(I1)/FN)  *CDNJG  (Z  (11) /FN)  ) 
COYZ  (I1)=CCYZ  (11)  +CABS  (  (YY(I1)  /FN)  *CON JG  (Z  (11) /FN)  ) 
SZX  (11)  =SZX(I1)  +((Z(I  1)  /FN)  *C0!UG  (XX  (11)  /FN)  ) 
COZX(I1)=COZX  (11)  ♦  CABS  (  (Z  (11)  /FN)  *CONJG  (XX  (11)  /FN)  ) 

3  0    CONTINUE 
70    CONTINUE 
C  THIS    BLOCK   AVERAGES     THE    DATA    POINTS    ADDED    IN    BLOCK    30 

C  ABOVE    3Y     »NR«.    AT    T3IS    POINT    THE    POWER    CURVES     ARE 

C  CONVERTED    INTC    ?Dr:Z3     SPECTRAL    DENSITIES. 

C  THE    COHERENCES    AND    STOKES    PARAMETERS    ARE 

C  COMPUTED     IN     ALL    THREE    PLANES. 

DO    33    13=1  ,N 

S2XY  (13)=  (REAL  (SXY  (13)  )  /FNP)  *4. 

S3XY  (13)=  (AIM  AG  (SXY(I  3)  )  /FNP»  *'4  . 

S2YZ  (13)  =  (  PEAL  (SYZ  (13  )  )  /FNR)  *'4. 

S3YZ  (13  )=  (AIM  AG  (SYZ  (I  3)  )  /FNP.)  *4  . 

52ZX(I3)  =  (PEAL  (SZX  (13  )  )  /FNR)  *!*. 

S3ZX  (13)  =  (  AIM  AG  (SZX  (I  3)  )  /FNR)  *'4  . 

COXY  (13)=  (COXY  (13)  )  /SQPT(ZX  (13)  *ZY  (13)  ) 

CCYZ.(I3)  =  (COYZ  (13)  )  /SQRT(ZZ(I3|  *ZY  (13)  ) 

COZX  (13)  =  (COZX  (13)  )  /3QRT(ZX(I3)  *ZZ  (13)  ) 

ZX  (I3)=ZX  (13)  *T/FNR 

ZY(I3)  =ZY  (13)  *T/FNR 

ZZ  (I3)=ZZ  (13)  *T/FNR 

S1XY(I3)=  (  (ZX  (13)  -ZY  (13))  /T) *2. 

S1YZ(I3)=  (  (ZY  (13)  -ZZ  (13))  /T)  *2. 

S1ZX(I3)=  (  (ZZ  (13)  -ZX  (13))  /T)  *2. 

SOXY  (13)=  (  (ZX  (13)  f  ZY  (13))  /T)  *2. 

SOYZ  (13  )=  (  (ZZ  (13)  +ZY  (I3)|  /T)  *2. 

S0ZX(I3)=  (  (ZX  (13)  *ZZ  (13))  /T)  *2. 
33    CONTINUE 
C  THE    NEXT     DO    LOOP    CONVERTS    THE    3-FISLD    POWER    SPECTRAL 

C  DENSITY    INTO     DECIBELS    REFERENCED    TO     1     NANOTESLA    **2 

C  PER    HERTZ. 
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DO    32    1=1 , N 

ZX  (I)  =10.  *ALOG10(ZX(I  )  ) 

ZY  (I)  =10.  *ALOG10(ZY(I)) 

ZZ  (I)  =10.  *ALOG10(ZZ(I  )) 
32    CONTINUE 
C  THE    NEXT    BLOCK    ELIMINATES    THE    FI?5T    10    DATA    POINTS    30 

C  THAT    THE    FREQUENCY    IS     .04     TO    13    HZ. 

C  S1,S2,    AND    S3    ARE    NORMALIZED    B?    SO. 

C  DEGREE    OF     POLARIZATION,    ORIENTATION    ANGLE    ALPHA, 

C  AND    THE    FLLIPTICIIY    ARE    ALSO    COMPUTED. 

DC    950    IU=  1,31  32 

FREQ(TU)  =  FREQ  (IU+10) 

FRQ2(IU)  =FPQ2  (IU+10) 

ZX  (IU)=ZX  (IU+10) 

ZY  (IU)  =ZY  (IU+10) 

ZZ  (IU)=ZZ  (IU+10) 

COXY  (IU)  =CCXY  (IU+1  0) 

S1XY  (IU)=S  1XY  (IU+1  0)  /SOXY  (IU+m 

S2XY(IU)=S2XY  (14+10)  /SOXY  (14+10) 

S3XY(IU)=S3XY  (IU+10)  /SOXY  (IU+10) 

PXY(IU)  =SQRT(51XY(I4)  **2+S2XY(I4)  **2+S3XY  (IU)  **2) 

DEXY(IU)=0.5*ARSIN  (S3  XY  (14) /PXY  (14)  ) 

SOXY  (IU)  =  1  0.  *ALOG  10(3  OXY  (IU+10)  ) 

COYZ(IU)  =COYZ  (IU+10) 

S1YZ  (IU)=S  1YZ  (IU+10)  / SO YZ (14+13) 

S2YZ(IU)=S2YZ  (IU+1  0)  /SO  YZ  (IU*  1  }  ) 

S3YZ  (IU)  =  S3YZ  (IU+10)  / SO YZ (14+10) 

PYZ  (IU)  =SQRT(S1YZ(I'4)  **2+S2YZ  (14)  **2+S3YZ  (IU )  **2) 

DEYZ(I4)=0.5*ARSIN  (3  3  YZ  (IU)  /?  YZ  (IU  )  ) 

SOYZ  (IU)  =  1  0.  *ALOG10(SOYZ(IU+10»  ) 

COZX(IU)  =COZX(IU+10) 

S1ZX(IU)=S  1ZX  (IU+10)  /S0ZX(IU+1D) 

S2ZX(IU)  =S2ZX  (IU+1  0)  /SOZX  (14+13) 

S3ZX(IU)=S3ZX  (IU+10)  /SOZX  (IU+10) 
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c 
c 
c 

c 
c 
c 
c 


PZX  (14)  =SQRT(S  1ZX(I4>  **2+S2ZX  (14)  **2+S3ZX(I4)**2) 
DEZX(I4)=0.5*ARSIN  (5  3  ZX  (I  4)  /PZX  (I  4)  ) 
S0ZX(I4)  =1  0.  *ALOS10(SOZX  (14+10)  ) 
AXI  (14)  =AT  AN(S2XY(I^>  /S1XY  (14)  |  /2  . 
AYZ  (14)  =ATAN(S2YZ(I'4|  /S  1 Y  Z  (14)  )  /2  . 
AZX  (14)  =AT  AN(S2ZX(I4)  /S1ZX(I4)  |  /2. 
950    CONTINUE 

VEBSATEC    PLOTS 

N?TS=10./DELTAF    +1  . 

'NPTS'     DETERMINES    NUMBER    OF    POINTS    NECESSARY 

FOR    THE     0   TO     10    HERTZ    RANGE    TO    BE    PLOTTED. 

FOR    THE     'ITS'     AND    »  RTB1     VALUES    REVIEW    THE    WRITE-UP 

FOR    THE     SUBROUTINE     PROCEDURE    • DEAWP' . 

ITB  (3)  =7 

TIB  (4)  =5 

ITB  (7)  =1 

ITB  (12)  =1 

RTB  (1)  =0.  0 

RTB  (2)=0.  0 

RTB  (3)  =ALAB(1) 

READ(5,3000)  TITLE 

CALL    DRAW?  (NPTS, FRQ2,  ZX, ITB, RT3) 

RTE  (3)=ALAB(2) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FRQ2r  ZY,  ITB,  RIB) 

RTB  (3)  =ALAB(3) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (NPTSrFRQ2,  ZZ,ITB,RTB) 

RTB(3)=ALAE(4) 

READ  (5,30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  CCXZ  , ITB,  3TB) 

RTB  (3)=ALA3(5) 
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HEAD  (5,  30  0  0)  TITLE 

CALL    DRAWP  (N?TS,FRQ2,  SOXf , ITB, RTB) 

RTB  (3)  =  ALAB(6) 

READ  (5,3000)  TITLE 

CALL    DRAWP  (NPTS,?RQ2r  S1  XY  ,ITB  ,  RTB) 

RTB  (3)=  ALA  B  (7) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  S2XY  , ITB, RTB) 

RTB  (3)  =ALAB(8) 

READ  (5,  3000)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  S3  XY  ,  IT3  ,  B  I  B) 

RTB  (3)=ALAE(4) 

READ  (5r  30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  COY2  rITB#  RTB) 

RTB  (3)  =ALAB(5) 

READ(5,3000)  TITLE 

CALL    DKAWP  (NPTS,?RQ2,  SOYZ  ,ITB,  RTB) 

RTB  (3)=  ALA  3(6) 

READ  (5,  30  0  0)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  S1YZ  ,  ITS,  RTB) 

RTB  (3)  =ALA3(7) 

READ  (5,30  0  0)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  S2YZ , ITB, R IB) 

RTB  (3)=ALAE(8) 

READ  (5,  3000)  TITLE 

CALL    DRAWP  (NPTS,FRQ2,  S3YZ  , ITB, RTB) 

RTB(3)  =ALAB(U) 

HEAD  (5,  3000)  TITLE 

CALL    DRAW?  (N?TS,FRQ2,  COZX  , ITB,  RTB) 

RT3(3)=ALAE(5) 

READ(5,3000)  TITLE 

CALL    DRAW?  (NPTS,FRQ2,  SOZX  ,  ITB, 3TB) 

RTB  (3)=ALAB(6) 

READ(5r3000)  TITLE 
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CALL    DRAWP  (NPTS,FRQ2, 
RTE  (3)=ALAB(7) 
READ  (5r  30  0  0)  TITLE 
CALL    DRAW?  (NPTS,FRQ2, 
PTB  (3)  =ALA3(8) 
READ  (5,  3000)  TITLE 
CALL    DRAW?  (N?TS,FRQ2, 
PTB(3)=ALAE(9) 
READ  (5,  30  0?)  TITLE 
CALL    DRAWP  (NPTS,FRQ2, 
RTE  (3)  =  ALA  3(10) 
READ(5,3000)  TITLE 
CALL    DRAW?  (N?TS,FRQ2, 
RTE  (3)=ALAE(1  1) 
READ(5f  3000)  TITLE 
CALL    DRAWP  (N?TS,FRQ2, 
RTB  (3)  =ALAB(9) 
READ  (5,30  00)  TITLE 
CALL    DRAW?  (NPTS,FRQ2, 
RTB  (3)=ALAE(10) 
READ  (5,  30  0  0)  TITLE 
CALL    DRAWP  (N?TS,FRQ2, 
RTB  (3)  =ALAB(1  1) 
READ(5,3000)  TITLE 
CALL    DRAWP  (N?TS,FRQ2, 
RTB  (3)=ALA3(9) 
READ(5f  30  0  0)  TITLE 
CALL    DRAWP  (N?TS,FRQ2, 
RTB(3)  =  ALAB(10) 
READ(5r3000)  TITLE 
CALL    DRAWP  (N?TS,FRQ2, 
RTE  (3)=ALAE(1  1) 
R3AD(5,  3000)  TITLE 
CALL    DRAWP  (NPTS,FRQ2, 


S1ZX, ITB, 3TB) 


S2ZX,ITB,RTB) 


S3ZX,ITB,RTB) 


AXY,  ITB,  RTB) 


?XY,IT3,RI3) 


DEX? rITB, SIB) 


AYZ,ITB,P.r3) 


PYZ,IT3,Rr3) 


DEYZ  , ITB, RTB) 


AZX,IT3r RIB) 


?ZX,IT3,Rr3) 


DEZX  ,  ITB, RTE) 
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3300    FORMAT  (6A8) 

STOP 

END 
//LKED.SYSUT1    DD    S? ACE= ( 1 3 2 4 ,  (200 , 233  )  ) 
//30.SYSIN    DD    * 

LA    MESA    SITS,     16    SEP    82,  1250-1330    LOCAL 

PSD    OF    X-COIL.     AMP    IN    D5  REF    NT**2/HZ    VS    HZ 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 

PSD    OF    Y-COIL.     AMP    IN    DB  RSF    NT**2/HZ    7  5    HZ 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 

PSD    OF    Z-COIL.    A KP    IN    DB  3.2  F    NT**2/HZ    VS    HZ 

LA    MESA    SITE,     16    SEP    32,  1250-1330    LDCAL 

COHERENCE    OF    X    B    Y    COILS  VS    FREQ  (HZ) 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
S0(10LOG)     VS    FREQ(HZ)  X-Y    PLANS 

LA    MESA    SITE,     16    SEP    32,  1250-1333    LDCAL 
S1/S0    VS    FREQ  (HZ)  X-Y    PLANE 

LA    MESA    SITE,     16    SEP    32,  1250-1330    LDCAL 
32/50    VS    FREQ  (HZ)  X-Y    PLANE 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
S3/S0    VS    FREQ  (HZ)  X-Y    PLANE 

LA    MESA    SITE,    16    SEP    82,  1250-1330    LDCAL 

COHERENCE    OF    Y    S    Z    COILS  VS     FREQ  (HZ) 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
S0(10LOG)     VS    FP.EQ(HZ)  Y-Z    PLANE 

LA    MESA    SITE,    16    SEP    82,  1250-1330    LDCAL 
S1/S0    VS    FREQ(HZ)  Y-Z    PLANE 

LA    MESA    SITE,    16    SEP    82,  1250-1330    LDCAL 
S2/S0    VS    FREQ  (HZ)  Y-Z    PLANE 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
S3/S0    VS    FREQ (HZ)  Y-Z    PLANE 

LA    MESA    SITE,    16    SEP    82,  1250-1330    LDSAL 

COHERENCE    OF    Z    5    X    COILS  VS     FREQ  (HZ) 

LA    MESA    SITE,    16    SEP    82,  1250-1330    LDCAL 
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S0(10LOG)     VS    FREQ(HZ)  Z-X    PLANE 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LOCAL 

S1/S0    VS    FREQ(HZ)  Z-X    PLANE 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LOCAL 

S2/S0    VS    FREQ(HZ)  Z-X    PLANE 

LA    MESA    SITE,    16    SEP    82,  1250-1330    LDCAL 

S3/S0    VS    FREQ(HZ)  Z-X    PLANS 

LA    MESA    SITS,     16    SEP    32,  1250-1330    LDCAL 
AXY    VS    FHEQ(HZ) 

LA    MESA    SITE,     16    SEP    82,  12  50-1330    LOCAL 
DEGREE    OF    POLARIZATION    VS    ?REQ(HZ)          X-Y 

LA    MESA    SITE,     16    SEP    32,  1250-1330    LDCAL 

ELLITICITY    VS    FREQ(HZ)  X-Y    PLANS 

LA    MESA    SITE,     16    SEP    32,  1250-1330    LDCAL 
AYZ     VS    FREQ(HZ) 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
DECREE   OF    POLARIZATION    VS    FREQ(HZ)          Y-Z 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LOCAL 

ELLITICITY    VS    FREQ(HZ)  Y-Z    PLAN'S 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
AZX    VS    FF.EQ(HZ) 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 
DECREE   OF    POLARIZATION    VS    FREQ(HZ)           Z-X 

LA    MESA    SITE,     16    SEP    82,  1250-1330    LDCAL 

ELLITICITY    VS    FREQ(HZ)  Z-X    PLANS 

/* 

//GO. FT20F001  DD    UNIT=3U03- 4, V3L=SSR=3MDT52, DIS?= (OL D, KEEP) , 

//  LABSL=  (1,NL,  ,  IN)  , 

//  DCE= (RECFM  =  F3,LRECL  =  3  2,  BLKSIZS=51 2 , DSN  =  2) 

/* 

// 
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APPENDIX    B 


CHEW'S    RID3E     Sir E    PREPARATION 


Chew's  Ridge  is  the  site  of.  previous  land-side 
geomagnetic  data  collection.  It  is  located  approximately 
twenty-three  air  miles  southeast  of  the  Naval  Postgraduate 
School  or.  a  hilltop  3,300  feet  abovs  sea  level.  It  was 
chosen  for  its  remoteness  from  the  local  power  grid  (3-4 
miles)  and  its  virtual  "line  of  sight"  RF  path  to  the 
school.  A  detailed  description  of  the  physical  layout  and 
operation   of   the   site   can   be   found   in    [Ref.    6]. 

The  original  intention  of  this  rssearch  was  to  collect 
data    from   this   remote   sits    to   provide      the    simultaneous   data 


necessary    equipment      continued    through    most    o: 


he  projec 


to    the   Naval 
rcvsd   impossible 


The   ca~a    transmission   from   this    s: 
Postgraduate  School  receiving  location 

with  the  transmitter/receivers  available. 

Previous  data  collected  from  Chew's  Ridge  was  collected 
using  a  Motorola  system  that  suffered  from  overheating 
failures  when  reguired  tD  operate  in  a  continuous  mode. 
Intially  another  attempt  was  made  to  repair  the  Motorola 
unit  but  after  one  last  failure  this  proved  economically 
impossible  and  if  was  felt  that  this  unit  did  not  have  the 
bandwidth  to  accomadate  the  addition  of  the  PCM  to  the 
system.  Consequently  a  RSPCO  :iodei  910  transmitter-receiver 
pair  was  ordered.  These  on  arrived  the  first  of  September 
and  initial  testing  showed  that,  although  they  were  capable 
of  operating  for  long  periods  of  time,  they  did  not  have  the 
required  20  KHz  bandwiith  to  accomadate  the  PCM  data. 
Efforts  are  currently  underway  to  expand  the  bandwidth  to 
meet  the  requirements. 
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Figure  B.1    Geographical  Area  Df  Data  Collscrion 
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To  properly  orient  tha  sensing  cdIIs  at  -this  site  it  is 
necessary  to  have  an  accurate  known  direction.  In  order  to 
accomplish  this,  inguiriss  were  mads  to  The  National  Fores- 
Service  to  determine  if  they  have  survey  data  for  the  area. 
Whan  it  proved  nonexistent,  a  survey  by  resection  was 
conducted.  A  theodolite  was  used  to  measure  relative  angles 
to  several  prominent  terrain  features  in  the  area  from  the 
Chew's  Ridge  Site.  These  measurements  were  made  three  times 
and    then      averaa=d.       The    iveraae      relative    ancrles      were   -her. 


drawn   o: 


siece   of    overlay    paper   and    oriented   on    a    1:50,000 


map  of  the  area.  This  allowed  for  the  best  fit  of  the  data 
and  prevented  blunders  in  the  recognition  of  the  terrain 
features.  The  most  prominent  terrain  feature  (Uncle  Sam 
Mountain)  was  then  chosen  for  the  end  of  the  orisnting  line. 
Tha  Grid  Azimuth  to  Uncle  Sam  Htn.  was  then  computed  from 
the  coordinates  of  the  two  locatiDns,  Finally  the  True 
Azimuth  was  calculated  to  be  257.4  degrees.  Estimated 
accuracy   of  this    measurement    is    one    degree. 
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